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Abstract: We compute the 3-loop 0{aas) correction to the Higgs boson production cross 
section arising from hght quarks using an effective theory approach. Our calculation probes 
the factorization of QCD and electroweak perturbative corrections to this process. We 
combine our results with the best current estimates for contributions from top and bottom 
quarks to derive an updated theoretical prediction for the Higgs boson production cross 
section in gluon fusion. With the use of the MSTW 2008 parton distribution functions that 
include the newest experimental data, our study results in cross sections approximately 
4 — 6% lower for intermediate Higgs boson masses than those used in recent Tevatron 
analyses that imposed a 95% confidence level exclusion limit of a Standard Model Higgs 
boson with Mh = 170 GeV. 
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1. Introduction 



The Higgs boson is the last undiscovered particle of the Standard Model. The search 
for the Higgs is a primary goal of the LHC program, and it is also a top priority at the 
Tevatron. A discovery of the Higgs boson is feasible with good confidence at the LHC for 
all mass values where the Standard Model remains consistent. The Tevatron experiments 
are becoming sensitive to Higgs signals in the H ^ WW channel for masses near the 
threshold Mh ~ 2M\y- Recently, the Tevatron collaborations reported a 95% confidence 
level exclusion of a Standard Model Higgs boson with a mass near Mh = 170 GeV [1,2]. 

Understanding the theoretical prediction is crucial to both the search for and exclu- 
sion of the Standard Model Higgs boson. Backgrounds to the Higgs signal are severe 
in many channels, particularly when a mass peak cannot be reconstructed such as in 
H — > WW — > Ivlv, and knowledge of the signal shape and normalization is needed to 
optimize experimental searches. Measurements of Higgs boson couplings will also require 
the best possible theoretical predictions [3,4]. The dominant production mode at both the 
Tevatron and the LHC, gluon fusion through top-quark loops, receives important QCD 
radiative corrections [5-7]. The inclusive result increases by a factor of 2 at the LHC and 
3.5 at the Tevatron when perturbative QCD effects through next-to-next-to-leading order 
(NNLO) are taken into account [8-10]. The theoretical uncertainty from effects beyond 
NNLO is estimated to be about ±10% by varying renormalization and factorization scales. 
A better perturbative convergence and a much smaller scale variation are found when typi- 
cal experimental cuts which suppress associated jet radiation at high transverse momentum 
and enhance the H — > WW signal at Tevatron and the LHC are implemented [11-16]. 

The importance and success in taming the QCD corrections to Higgs production have 
shifted attention to electroweak corrections to the Higgs signal. The authors of Refs. [17,18] 
pointed out important 2-loop light-quark effects; these are pictured in Fig. (1) of this 
manuscript and involve the Higgs coupling to W- or Z-bosons which then couple to glu- 
ons through a light-quark loop. These terms are not suppressed by light-quark Yukawa 
couplings, and receive a multiplicity enhancement from summing over the quarks. A care- 
ful study of the full 2-loop electroweak effects was performed in Ref. [19]. They increase 
the leading-order cross section by up to 5 — 6% for relevant Higgs masses. However, it is 
unclear whether these contributions receive the same large QCD enhancement as the top 
quark loops. If they do, then the full NNLO QCD result is shifted by -|-5 — 6% from these 
electroweak corrections. If not, this 5 — 6% increase from light quarks would be reduced 
to 1 — 2% of the NNLO result. As this effect on the central value of the production cross 
section and therefore on the exclusion limits and future measurements is non-negligible, it 
is important to quantify it. The exact computation of the mixed electroweak/QCD effects 
needed to do so requires 3-loop diagrams with many kinematic scales, and 2-loop diagrams 
with four external legs for the real-radiation terms. Such a computation is prohibitively 
difficult with current computational techniques. 

In this paper we compute the QCD correction to the light-quark terms in the Higgs 
production cross section using an effective theory approach. We justify our approach 
rigorously by applying a hard-mass expansion procedure to the full 3-loop corrections. 
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This technique reduces the calculation to the evaluation of 3- loop vacuum bubbles. The 
effective theory is formally valid only for Mh < M\y- However, there are reasons to believe 
that the if-factor computed with the effective theory has an extended range of validity. In 
the top-quark contribution to gluon fusion, the effective theory obtained after decoupling 
the top quark is formally valid only for Higgs boson masses Mh < 2mt. Nevertheless, the 
ET-factor obtained is an extremely good approximation to the exact one for Higgs boson 
masses up to Mh ~ 1 TeV [5] . We find that the correction to the light-quark terms is not 
as large as those affecting the top-quark contribution. Nevertheless, the two corrections 
have the same sign, and the numerical effect of the difference is small, indicating that the 
5 — 6% shift is indeed realized. 

A second goal of this manuscript is to provide the most up-to-date QCD prediction 
for the Higgs boson production cross section in gluon fusion for use in setting Tevatron 
exclusion limits. The CDF and DO collaborations [20] currently use results from Ref. [21], 
which are several years old, augmented by the light-quark corrections from Ref. [18]. The 
following aspects of the analysis given must be updated to account for recent developments. 

• Both CTEQ and MRST parton distribution functions (PDFs) have shifted signifi- 
cantly in the past several years due to an improved treatment of heavy-quark mass 
effects at low [22-25], and inclusion of several Tevatron Run 11 data sets [24]. 
The new PDFs have a different as{Mz) and gluon distribution, and decrease the 
predicted production cross section. 

• The analysis in Ref. [21] used the if-factor computed in the effective theory with 
the top quark integrated out for both the top- and bottom-quark contributions. The 
NLO QCD correction to the bottom-quark contribution is known to be smaller than 
the NLO top-quark i^-factor [5,26]. This effect increases the predicted production 
cross section. 

• An updated treatment of the 2-loop light-quark contributions from Ref. [19], together 
with the QCD correction to these terms evaluated here, leads to a slightly smaller 
increase near the Mh ~ 2Mw threshold than used in the Tevatron analysis. 

We present results for the Higgs boson cross section accounting for these effects. We 
account for the effect of soft-gluon resummation at the Tevatron by presenting values for 
the scale choice hf = IJ-r = Mh/^,, which is known to very accurately reproduce the 
reference value of the resummation result [21] for a wide range of Higgs boson masses, 
and provide an estimate of the remaining theoretical uncertainties arising from unknown 
higher-order terms and PDF errors. The updated numerical values for the cross section are 
approximately 4 — 6% lower than those used in Tevatron analyses for Higgs boson masses 
in the range 150 — 180 GeV, and motivate a reanalysis of the Tevatron exclusion limits. We 
present a detailed discussion of the uncertainties arising from scale variation, PDF errors, 
and other theoretical effects. 

Our paper is organized as follows. In Section 2 we describe our calculation of the 3- 
loop light-quark correction to the Higgs production cross section, detailing the formulation 
of the effective theory and technical aspects. In Section 3 we present numerics for both 
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the light-quark electroweak shifts and the updated inclusive cross section. We conclude in 
Section 4. 



2. Calculational Details 



The cross section for Higgs boson production in hadronic collisions can be written as 

cr(s,M^) = V/ dxi [ dx2fi/hA^i,IJ'F)fj/h2ix2,IJ'F) [ dzdfz-^^^ 
■ ■ Jo Jo Jo \ X1X2S/ 



zuij {z; as{iJ.ji),aEW,M^/iij^;M^/iij,) 



(2.1) 



Here, y/s is the center-of-mass energy of the hadronic collision, fiji and respectively 
denote the renormalization and factorization scales, and the fi/^ denote the parton den- 
sities. The quantity za is the partonic cross section for the process ij ^ H + X with 
h 3 = 5) 9; Q- As indicated, it admits a joint perturbative expansion in the strong and 
electroweak couplings. 

The leading term in the partonic cross section arises from a one-loop correction coupling 
the Higgs boson to gluons via a top-quark loop. Considering only QCD interactions for the 
moment and suppressing the scale dependences for simplicity, the partonic cross section 
can be written as 

ai,=(7(0)Gi,(z;a,), (2.2) 



with 



a 



(0) 



QH = myMfj, and 



512V27r 

2-{l-AqH)H[-r, 

1, 



-r; 

QH 



-r;x) = ^In^ 



\/4 + X - ^/x 



(2.3) 
(2.4) 

(2.5) 



V4 + X + 

The contribution from bottom quarks in the Standard Model is also numerically relevant; 
we discuss its inclusion later in this paper. The coefficient functions can be expanded in 
the strong coupling constant as as 



00 

G,j{z;as) = J2{^yG^\z), 

n=0 



(2.6) 



with the leading term given by Gf^\z) = 5ig6jg5{l — z). The NLO coefficient functions have 



been computed in Ref. [5] retaining the exact dependence on the quark mass. The NNLO 
results in the large ruq limit, relevant when 2mq > Mh, were derived in Refs. [8-10]. The 
NLO result in this limit was obtained in Refs. [6,7]. The QCD corrections have a large 
effect on the predicted cross section, increasing it roughly by a factor of 2 at the LHC and 
by a factor of 3.5 at the Tevatron. 
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Important electroweak corrections arise from two-loop diagrams containing an internal 
quark loop where the Higgs boson couples to W- and Z-bosons. An example diagram is 
shown in Fig. (1); we henceforth refer to these corrections as light-quark electroweak con- 
tributions, while the quark Yukawa coupling dependent terms discussed above are denoted 
as heavy-quark QCD contributions. The light-quark diagrams are not suppressed by quark 



Yukawa couplings, and therefore have a multiplicity enhancement from summing over light 

-,(0) 



quarks. The inclusion of these contributions modifies the term proportional to g[^\z) in 



Eq. (2.2). The partonic cross section becomes 



n=l 

with 
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Qu + Qt 



(2.8) 



Qfj is the expression for the two-loop light-quark contributions; its analytic form in terms 
of generalized harmonic polylogarithms can be found in Ref. [17]. A calculation of the 
corrections with the light-quark loop replaced by a top-quark, or the top-bottom doublet 
in the case of the M^-boson, was first performed in Ref. [27]. A careful numerical study of 
these electroweak corrections utilizing the complex-mass scheme to handle the threshold 
regions Mh 2 Mw^z was performed recently in Ref. [19]; this study also includes effects 
from internal top quarks coupling to the W and Z. The full corrections increase the 
leading-order cross section by +5 — 6% for Higgs boson masses in the range 120 — 160 GeV. 



Figure 1: Example two-loop light-quark diagram contributing to the Higgs boson production cross 
section via gluon fusion. 




The cross section in Eq. (2.7) includes corrections to the leading-order result valid 
through 0{a) in the electroweak couplings and to 0{a'^) in the QCD coupling constant in 
the large top- mass limit upon inclusion of the known results for . Since the perturba- 
tive corrections to the leading-order result are large, it is important to quantify the effect 
of the QCD corrections on the light-quark electroweak contributions. This would require 
knowledge of the mixed 0{aas) corrections, which arise from 3-loop diagrams. In lieu of 
such a calculation, the authors of Ref. [19] studied two assumptions for the effect of QCD 
corrections on the 2-loop light-quark diagrams. 

• Partial factorization: no QCD corrections to the light-quark electroweak diagrams 
are included, so that the cross section is given by the expression in Eq. (2.7). With 
this assumption, electroweak diagrams contribute only a +1 — 2% increase to the 
Higgs boson production cross section. 
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• Complete factorization: the QCD corrections to the electroweak contributions are 
assumed to be identical to those affecting the heavy-quark diagrams, and the partonic 
cross section is therefore taken to be 

a?/ = c7^EwGij{z-as) (2.9) 

with the full QCD coefficient function multiplying both the heavy- and light-quark 
contributions. In this case the light-quark diagrams increase the full NNLO QCD 
production cross section by -|-5 — 6%. 

The resulting shift in the central value for the Higgs boson production cross section can 
have a non-negligible effect on exclusion limits at the Tevatron, and on future exclusion 
limits or measurements of Higgs boson properties at the LHC. 

Wc discuss later in this manuscript the inchision of bottom-quark contributions to the 
Higgs production cross section. We define for future reference the following point-like cross 
sections: 



(0) _ Gpa 



2 



gj^[2Re(fiftV)]. (2.10) 



denotes the squared bottom-quark loop, af*^ the interference between the top and 

bottom loops, and af'^j: the interference between the top-quark contribution and the light- 
quark terms. 

2.1 The Effective Lagrangian formulation 

A rigorous test of factorization of QCD and electroweak corrections to Higgs boson pro- 
duction in the Standard Model for all values of Mh would require a full 3-loop calculation 
containing several mass scales. Such a computation is seemingly beyond current technical 
capabilities. However, in the limit Mh < Mw,z, an approximate result can be obtained 
by expanding around the point Mh/Mw,z = 0. This is the same approach used to derive 
the heavy-quark result when Mh /{2mt) < 1. Although experimentally Mh > M\y_z, it 
is known that the approximate NLO correction to the heavy-quark result matches the ex- 
act NLO correction extremely well up to Mh ~ 1 TeV due to the structure of the QCD 
corrections. This provides some reason to believe the same is true for the light-quark 
contributions. 

The Mh = expansion is most clearly formulated using an effective Lagrangian, 
and we review this approach here. The leading term in the expansion for the top-quark 
contribution in Eqs. (2.3) and (2.4) can be obtained via Feynman rules derived from 

C^ff = -a.^HGl^G'^'^r (2.11) 
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The Wilson coefficient arising from integrating out the heavy quark is [28, 29] 
Ci = (1 + a,Ci, + a'iC2q + 0(a^)} , 

where Ug = ag/ir, Np = 5 is the number of active quark flavors and Lt = ln(;u|j/mj). We 
now include the leading term in the Mh /Mw^z expansion in the Wilson coefficient. It can 
be obtained by expanding Qfj^ given in Ref. [17]; denoting the contribution by Xew, we 
find 

Ci = (1 + Xew + asCiq + a^Ca^ + 0(4)} , 



Xew 



3a 



5 '^2 4 

4 - 3^iy + 



+ 4 , (2.13) 



where sw, cw are respectively the sine and cosine of the weak-mixing angle. 

The QCD corrections modify the Wilson coefficient to include terms of 0{Xew^s) and 
C(AEvi/a^). The extent to which factorization of electroweak and QCD corrections holds 
becomes a question regarding to what extent the Wilson coefficient can be written as a 
product of separate QCD and electroweak factors. We denote the exact coefficients of these 
terms as Ciw and C2w respectively, and introduce below in Eq. (2.14) the exact Wilson 
coefficient and the f actor ized hypothesis: 

Cl = {l + ^EW [l + asCiw + a^C2«,] + QgCiq + a1C2q} , 

Ct' = (1 + >^Ew) {1 + o.Ci, + alC2g} . (2.14) 

Factorization holds if Ci^j = Ciq and C2w = C2q- We will derive here the Ciw coefficient 
by expanding the 3-loop QCD corrections to the light-quark electroweak diagrams to test 
this. We do not compute C2W1 but will study later the numerical effect of various choices 
for this term. 

2.2 Calculational approach 

We begin by generating all 3-loop diagrams for g(pi) + g{p2) H(ph) containing two 
internal W- or Z-boson propagators coupling to a light-quark loop. Examples are shown 
in Fig. 2. Contributions with internal top quarks attached to the W 01 Z are very small 
at the 2-loop level for the Higgs boson masses we consider, as can be seen from Ref. [19], 
and can be safely neglected. For the Z-boson and for a single quark flavor, there are 51 
such non- vanishing diagrams. Examples are shown in Fig. (2). The only difference for 
the VF-boson is the change of quark flavor at the vertex. All such details are accounted 
for in Xew Siud do not affect the computation of Ciw We perform a Taylor expansion 
of the integrand of each diagram in the external momenta pi,2- This is most conveniently 
performed by applying the following differential operator to each diagram [30] : 

00 

n=0 
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where T denotes a diagram. The first few terms in the sum are 

i^o = l, i^i = ^Di2, ^2 = -2(^3^^^{DnD22-dDL}, (2.16) 

where = g^-^- The amphtude arising from summing over Feynman diagrams can 
be written as 

^ J- = ^ L^, - e^M^UP2) = Mf!le^M^UP2)- (2.17) 

L Pl-P2 ) 

The coefficient .4 can be obtained by acting with the appropriate projection operator: 

The leading term in the expansion of A must be finite, and gives Ciw upon comparison 
with Eqs. (2.11) and (2.14). Sub-leading terms in the Taylor expansion need not be finite. 
We note that the 2-loop light-quark contribution has a cut first at = M"^ ,^ because of 
helicity fiow along the internal quark line [27] , indicating that the radius of convergence of 
the expansion is Mh < M\y,z- The convergence is unchanged when the 3-loop corrections 
are added. 

The validity of the effective theory and the Taylor expansion we utilize for Mh < Mw,z 
is most clearly seen by applying the hard-mass expansion procedure [31] to the 3-loop 
diagrams. This technique allows one to derive the asymptotic limit of a diagram in the 
limit of large internal masses via an expansion in subgraphs: 

Jc-r^^^r/7°^fc,P,-^7- (2.19) 

7 

In this expression, V denotes all loop-momenta dependent pieces of the diagram JT. 7 
denotes the subgraphs, which are defined as those combinations of internal lines that con- 
tain all propagators with the heavy mass M^^^ and that are one-particle irreducible with 
respect to the massless lines. T^^p. indicates the Taylor expansion of the subgraph with 
respect to the external momenta pi and also the loop momenta k that are external to the 
subgraph. is ^ reduced graph, which is what remains of the diagram upon removal of 
a subgraph. 

To prove the validity of our procedure, we must show that the only subgraphs con- 
tributing to the leading Myy^z term are the full 3-loop diagrams themselves, and the 2-loop 
diagrams of Fig. 1 that give Xew- If this occurs, then all contributions are produced by 
the effective Lagrangian of Eq. (2.11). The 3-loop terms give Ci^,, while the 1-loop reduced 
graphs multiplying the 2-loop subgraphs are given exactly by the first quantum corrections 
in the effective theory. This is indeed what occurs. We sketch briefly below the steps of 
the calculation. 

On dimensional grounds the leading term of the coefficient A appearing in Eq. (2.17), 
which comes from summing all 3-loop graphs, must scale as ^ ~ gnvv^H/^w Zi '^^ere 
9HVV is the HVV coupling that has mass dimension one. The subgraph obtained by 
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Figure 2: Example three-loop light-quark diagrams contributing to the Ci^ term in the Wilson coef- 
ficient. 



expanding only the two massive gauge boson propagators goes as qhw l^w ™^ does 
not contribute to the leading term. The possible 1-loop subgraphs contain either the Higgs 
boson coupling to two quarks, or two quarks and a gluon. These subgraphs contribute only 
at Qhw l^w z when summed. We similarly find that the only 2-loop subgraphs are those 
contained in Xew-, which shows the validity of our approach. 

In addition to the explicit check outlined above, we make two more remarks that indi- 
cate the validity of the effective theory. The divergences of the effective theory match the 
universal structure of infrared divergences as given by the Catani factorization formula [32] . 
Also, the only other operator in the effective Lagrangian describing Higgs interactions with 
gluons and massless quarks that could contribute at order gHW /^w z Hqpq [29]. 
When inserted into a loop of light quarks and coupled to two gluons, this contribution 
gives scaleless integrals and vanishes; we note that this is a simple confirmation of our 
statement above that the only 1-loop reduced graphs that contribute are those generated 
by the effective Lagrangian of Eq. (2.11). 

We now proceed with our calculation. All integrals appearing in A can be mapped to 
the following topology: 



1 



A-2 - /l/f2 
'^3 ^^■^W,Z 



^3 



[(fci - k2Yr [{k2 - k:,Yr [(^3 - kiY. 



V6 ■ 



(2.20) 



These can be reduced to a small set of master integrals using what is by now standard 
technology based on the solution of integration-by-parts identities [33-35]. We find that 
only the integrals T(l, 0, 1, 1, 1, 0) and T(l, 1, 1,0, 1, 1) are needed to describe this process. 
It is straightforward to express both as a simple product of Gamma functions. 

After a computation following the approach outlined above, we obtain the primary 
analytic result of this paper: 

Ci^ = I- (2.21) 
b 

Two points should be noted regarding the comparison of this with the factorization hy- 
pothesis Cl^" = Cig = 11/4. First, there is a fairly large violation of the factorization 
result: (Cig — Ciw)/Ciw ~ 1-4. However, both expressions have the same sign, and a large 
difference from the +5 — 6% shift found before does not occur. We study in the next section 
the numerical effect of Ci^ and the unknown C2w 



-8- 



3. Numerical Results 



We present here numerical results for the Higgs production cross section and the shift 
arising from light-quark clcctroweak diagrams including the effect of Ciw and C2w Our 
purpose is two- fold: to study the numerical effect of the correction computed in the previous 
section, and to provide an updated prediction for the inclusive cross section for use in 
experimental studies. We include all currently computed perturbative effects on the cross 
section. These include the NNLO K-factor computed in the large-mt limit and normalized 
to the exact m^-dependent LO result, the full light-quark electroweak correction and the 
0{as) correction to this encoded in Ci^, and the bottom-quark contribution. We define 
the following cross sections for use in our discussion: 



^QCD — ''^ij {Z,as) + (Tjj U-j {Z) + a^ f^U^j [Z)Kti, , 

JO) ^(0) 



^LO _ JO) ^(0) . X 



NLO 
'^EW 



+ 0-s{Clw — Ciq) + a1 {C2w — C2q + Ciq{Ciq — Ciw)] 

+asG[f {z) [1 + a,(Ci^ - Ci,)] + a^G^f } , 

NNLOCF _ (0) fi.r^.^ a 

„best NNLO . NNLO /o -i\ 

^ = (^QCD + (^EW ■ [■^■^) 

Unless noted otherwise, all results use the MSTW 2008 distribution functions [24] at the 
appropriate order noted in the superscript of the cross section. We briefly describe here 
the content of these several terms, (^qcd^ includes contributions from both top- and 
bottom-quark loops, with cj^'^'' and af^ defined in Eq. (2.10). The QCD corrections to 
the top-quark in the large-mj limit are encoded in {z; a^). The NLO iC-factors for the 
squared bottom-quark term and the interference term as derived from Ref. [26] 
are included in Kbi, and Kth respectively. We note that both the study in Ref. [21] and the 
Tevatron analysis put the bottom-quark terms in (7^^\ and therefore use the same iT-factor 
for both top- and bottom-quark loops. This results in an underestimate of the cross section, 
since the effect of these terms is negative; while the NNLO i^— factor for the top-quark 
term is roughly 2.1 and the NLO iC-factor is roughly 1.8 with MSTW2008 PDFs, iCbb and 
Kfi, only vary between 1.2 and 1.5 for Higgs boson masses between 120 — 180 GeV. 

The remaining terms in Eq. (3.1) denote the contributions from light-quark diagrams 
in various approximations, af'^j denotes the interference between the dominant top-quark 
term and the light quarks defined in Eq. (2.10); in our numerics we use the exact results of 
Ref. [19] which are valid for arbitrary Higgs boson masses. The squared light-quark term is 
numerically irrelevant, cr^^ includes only the 2-loop diagrams computed in Refs. [17, 19] 
and is equivalent to the partial factorization hypothesis defined in Sec. 2. cr^yp includes 
the 0{as) correction to these diagrams computed in the effective theory and encoded in 
Ciw o'ew^^ includes the full 0{al) correction to the light-quark diagrams including the 
unknown coefficient C2w We study numerically below various choices for C2w <^ew^'~' '^^ 
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is the complete factorization hypothesis defined in Sec. 2. Finally, a '' is the current 
best prediction for the Higgs boson cross section including all effects of top, bottom, and 
light quarks with the best estimates of their associated QCD corrections. To approximately 
implement the effects of soft-resummation, we make the scale choice fiR = fip = = Mh/2. 
This choice is known to reproduce the central value of the resummation results to better 
than 1% accuracy [21]. It has been pointed out that the choice /x = Mh may not correctly 
describe the Higgs production process [9, 12], and that the perturbative convergence is 
improved for // = mij/2 [36]. The effect of resummation is also smaller for this scale 
choice [37]. We evaluate the electroweak corrections using Gp, Mw and Mz as input 
parameters. We use the pole mass mt = 170.9 GeV for the top quark and the MS mass 
'nT-bit^R) for the 6-quark with the input value mb(10 GeV) = 3.609 GeV [38] . The choice 
of pole or MS mass for the top quark has a negligible effect on the numerical results. 

We comment here on the numerical validity of the large-nit limit for the top-quark 
squared contribution. Factoring out the exact top-quark dependence and computing the K- 
factors in the effective theory, as we do here, gives an exceptionally accurate approximation 
to the full result. We have checked with an exact calculation at NLO that our result agrees 
to better than 1% in the kinematic range relevant for the Tcvatron studies; confirmations of 
this result have been obtained using several independent codes [39]. Finite mt corrections 
to the NNLO coefficient function have been shown to affect the i^-factor by less than 1% 
when the full mt dependence is factored out [40,41]. We note that we compute exactly 
at NLO the top-bottom interference and bottom-bottom diagrams. We conclude that the 
error arising from our treatment of the top quark mass is at the percent level or less. 

We begin by studying the percentage shifts arising from electroweak effects on the 
Higgs boson production cross section at both the Tcvatron and the LHC in Fig. (3). The 
results shown in these plots are = 100 x (^ew/'^qcd^^ with the cross sections defined 
in Eq. (3.1). The close agreement between S^w^^ ^'^'^ ^ew"^'~^^ occurs because the 
differences C\w — C\q and C2w — in Eq. (3.1) are small compared to the effects of 
in (T^y^r^^ ■ Furthermore, the unknown C2w coefficient does not significantly alter 
the size of the electroweak shift. Our calculation confirms that the Higgs boson production 
cross section receives almost the entire numerical correction indicated by the complete 
factorization hypothesis. 

We now combine all effects into a best prediction for the Higgs boson production cross 
section, cr^^** defined in Eq. (3.1). We set C2w = to produce these numbers. As discussed 
previously, several updates must be performed to the analysis in Ref. [21] and therefore 
the Tevatron exclusion limits. 

• The -fC-factors for aj^^ and a^^jj are now known to be significantly smaller than those 
for top-quark term a^^^ [5,26]. For example, in the pole-mass scheme the NLO K- 
f actor for the top-bottom interference is approximately 1.3 for Mh = 150 GeV, as 
compared to over 1.7 at NLO and 3 at NNLO for the top-quark term; the K-factor 
for the bottom-quark piece is roughly 1.5 for this Higgs mass. The study in Ref. [21] 
utilized the top-quark X-factor for all three terms. As mentioned above we use 
the MS 6-quark mass instead of the pole mass of Ref. [21]. Our final cross section 
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Tevatron 



LHC 




Figure 3: Relative shifts to the Higgs boson production cross section at the Tevatron (left panel) and 
LHC (right panel) arising from light-quark diagrams. All curves are normalized to the full NNLO top- 
quark cross section and are produced for renormalization and factorization scales ^ir = fip = Mh/2. 
The various lines are described in detail in the accompanying text. 



numbers with the MS 6-mass are about 1.5% larger than the numbers computed 
in the pole scheme. Although the results in both schemes are very similar after 
the inclusion of NLO corrections, previous studies have shown some preference for a 
running mass [5,42] so we present numbers for this scheme. 

• Updated PDF extractions by CTEQ and MRST with an improved treatment of 
heavy-quark effects at low have a significant effect on the cross section. 

• The electroweak terms derived above must be added. 

These corrections all have an important effect on the predicted cross section. For illus- 
tration, we show below the sequential effect of making these changes on the cross section 
for Mh = 170 GeV. We begin by reproducing the fip = fiR = Mh/^ numbers given in 
Ref. [21], which also matchs the reference value for the resummed result, by implementing 
MRST 2002 NNLO PDFs, using the masses = 176 GeV and rub = 4.75 GeV, mul- 
tiplying all top- and bottom-quark terms by the ET-factor appropriate for the top quark, 
and removing the 2-loop light-quark terms. We then perform the following changes: (1) 
we switch to MRST 2008 NNLO PDFs; (2) we switch to the current extracted top-quark 
mass mt = 170.9 GeV, and to fhh, and use the NLO Ktb and Khb respectively to model the 
QCD corrections to the top-bottom interference and the bottom-quark squared contribu- 
tion; (3) we implement the electroweak corrections described above. The results are shown 
in Table 1. The effect of moving to MSTW 2008 PDFs is large and negative, although the 
other effects compensate to a large degree. We apply all of these corrections to provide 
up-to-date values for the inclusive Higgs boson production cross section in Table 2 for the 
scale choice // = -/Vf/f/2 that accurately reproduces the effect of soft-gluon resummation. 
These numerical values are 4 — 6% lower than values previously used by the Tevatron 
collaborations to establish exclusion limits on the Standard Model Higgs boson [20] for 
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original 


MSTW 2008 PDFs 


Ktb, Kbh 


EW effects 


0.3542 


0.3212 


0.3377 


0.3444 



Table 1: Shifts in the Higgs boson production cross section resulting from the changes detailed in the 
text. All numbers are In picobarns. 
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i.41 ( \± 1 /o pCtI j 


iOU 


n A'iAA /'-UQCJ?^ T-.rlf^ 
U.4o44 \±\) /Q PdlJ 


115 


1.243 (±7% pdf) 


165 


0.3854 (±9% pdf) 


120 


1.094 (±7% pdf) 


170 


0.3444 (±10% pdf) 


125 


0.9669 (±7% pdf) 


175 


0.3097 (±10% pdf) 


130 


0.8570 (±8% pdf) 


180 


0.2788 (±10% pdf) 


135 


0.7620 (±8% pdf) 


185 


0.2510 (±10% pdf) 


140 


0.6794 (±8% pdf) 


190 


0.2266 (±11% pdf) 


145 


0.6073 (±8% pdf) 


195 


0.2057 (±11% pdf) 


150 


0.5439 (±9% pdf) 


200 


0.1874 (±11% pdf) 


155 


0.4876 (±9% pdf) 







Table 2: Higgs production cross section (MSTW08) for Higgs mass values relevant for Tevatron, with 
A* = M-R = = The total cross section o-'^'*' is defined in Eq. (3.1). The theoretical errors 

PDFs are shown in the Table; the scale variation is roughly constant as a function of Higgs boson 

mass. Other potential sources of theoretical error are discussed in the text. 

Mh = 150 — 170 GeV, and motivate a reconsideration of their results. The cross section 
for the exclusion point is reduced by 6% from what was used in the Tevatron analysis. 

We also estimate the current theoretical uncertainties arising from uncalculated higher- 
order terms and PDF errors. To estimate the errors from higher-order effects we vary 
the scale fi in the range [M///4, M/f], which is a factor of two around the central value 
l^R = l^F = M = ^///2. For the PDF errors we use the error eigenvectors provided with the 
MRST 2008 fit. The scale errors are constant with Higgs mass to very good approximation, 
and are [—11%, +7%]. The PDF uncertainties have a slight dependence on the Higgs boson 
mass, as shown in Table 2. 

4. Conclusions 

In this paper, we considered mixed QCD-electroweak corrections to the Higgs boson pro- 
duction cross section in the gluon-fusion channel. Working in an effective field theory valid 
for mn < Mw, we provided the leading term of the three-loop 0{aas) contributions due 
to diagrams containing light quarks. This result allows us to check the factorization of elec- 
troweak and QCD corrections proposed in Ref. [18,19]. We showed that, despite a fairly 
large violation of the factorization hypothesis, a significant numerical difference from the 
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prediction of this hypothesis is not observed due to the structure of the QCD corrections. 
We combined the 2-loop hght-quark diagrams based on the complex-mass scheme for the 
W- and 2^-bosons [19] with our new 3- loop 0{aas) correction and with contributions from 
top and bottom quarks to provide an updated theoretical prediction for the production 
cross section. We found values that are 4 — 6% lower than those currently used by the 
Tevatron collaborations in the analysis that led to the 95% confidence level exclusion on a 
Standard Model Higgs boson with Mh = 170 GeV. Our results motivate a reconsideration 
of the Tevatron exclusion limits. 
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Note added: After this manuscript was first submitted, the MSTW collaboration released 
an update of their PDF distribution which includes new Tevatron Run II data, and is meant 
to supersede their previous releases [24]. The new distribution has a lower value of as{Mz) 
and uses the Tevatron Run II dijet data. The effect on the Higgs cross section of changing 
from MRST 2006 to MSTW 2008 is large; the cross section shifts downwards by nearly 
15% from MRST 2006. For illustrative purposes, we include below in Table 3 the effect 
on the total Higgs cross section for Mh = 170 GeV coming from the change in PDF sets 
from 2001 to 2008. The PDF uncertainties as estimated by the error eigenvectors have also 
increased significantly in the MSTW 2008 distribution. 



MRST 2001 


MRST 2004 


MRST 2006 


MSTW 2008 


0.3833 


0.3988 


0.3943 


0.3444 



Table 3: The Higgs production cross section in picobarns for Mh = 170 GeV at the Tevatron, using 
several different PDF distributions. 
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